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Abstract 

This paper characterizes when a Delone set X in W 1 is an ideal crystal in terms of restrictions 
on the number of its local patches of a given size or on the hetereogeneity of their distribution. 
For a Delone set X, let Nx(T) count the number of translation-inequivalent patches of radius 
T in X and let Mx(T) be the minimum radius such that every closed ball of radius Mx{T) 
contains the center of a patch of every one of these kinds. We show that for each of these 
functions there is a "gap in the spectrum" of possible growth rates between being bounded and 
having linear growth, and that having linear growth is equivalent to X being an ideal crystal. 

Explicitly, for Nx(T), if R is the covering radius of X then either Nx(T) is bounded or 
Nx(T) > T/2R for all T > 0. The constant 1/2R in this bound is best possible in all dimensions. 

For M X (T), either M X {T) is bounded or M X (T) > T/3 for all T > 0. Examples show that 
the constant 1/3 in this bound cannot be replaced by any number exceeding 1/2. We also show 
that every aperiodic Delone set X has Mx(T) > c(n)T for all T > 0, for a certain constant 
c(n) which depends on the dimension n of X and is > 1/3 when n > 1. 
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1. Introduction 

The discovery in 1984 of quasicrstalline materials, which are strongly ordered aperiodic 
structures, has generated renewed interest in precisely delineating the boundary between dis- 



crete sets that are fully periodic (ideal crystals) and those with less global order, see Radin [24] 
and several of the articles in , edited by Moody. This paper studies such questions in discrete 
geometry, investigating how strong the restrictions on local features of discrete point sets in W 1 
must be to enforce crystallinity. Such questions originally arose in geometric crystallography, 
see pf. The questions we study also have connections with the ergodic theory of M n -actions 



[17|; there are parallel questions in the combinatorics of words, in one or several dimensions, 



which we briefly consider. 

We study the patch-counting function Nx(T) and the repetitivity function Mx(T) of Delone 
sets X, where patches are identified only up to translation equivalence. (These concepts were 
studied in jl7| and are defined in §2.) We show that both functions have a "gap in the spectrum" 
of their possible growth rates: either they are bounded, which happens precisely when X is an 
ideal crystal, or they grow at least linearly in the radius T of the patch. Growth rates like log T or 
s/T are impossible for these functions. These results can be viewed as characterizations of ideal 
crystals: when X is an ideal crystal both functions are bounded, but nevertheless sufficiently 
slow linear growth of either of these functions is enough to characterize ideal crystals. They 
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can alternatively be viewed as giving linear lower bounds on the growth rates of these functions 
for Delone sets X that are not ideal crystals. 

The problem of characterizing ideal crystals in terms of restrictions on allowed types of local 
patches was first studied for "regular point systems" in 1976 by Delone et al. B for isometry- 
equivalence classes of local patches. They identified a radius T (depending on the dimension 
and the Delone constants of X) such that every Delone set X whose patches of radius T are 
all isometry-equivalent is a regular point system and therefore an ideal crystal. Dolbilin et al. 
0] gave an extension of this result to "multiregular point systems" , which comprise all ideal 
crystals: for each k > 1 they gave a radius T (depending on the dimension, the Delone constants 
of X and on k) such that if there are at most k isometry-equivalence types of patches of this 
radius, then X is an ideal crystal with at most k isometry-equivalence types of patches of any 
radius. This gives an effective way of deciding, for a given k, whether X has degree of regularity 
< k. It does not, however, give a way of deciding whether a given set X has some finite degree 
of regularity; that is, whether X is an ideal crystal. This latter drawback is unavoidable, 
because with any finite sample of X there is no way to tell that it is not part of an ideal crystal 
whose unit cell is larger than the size of the sample. However, the dependence of T on k in Q 
is linear, and it follows that there is a constant c (depending only on the dimension and the 
Delone constants) such that if there is some T with at most cT isometry patches of radius T 
then X is an ideal crystal. 



For the patch-counting function Nx(T), our result (Theorem 2J.) is an analogue of this 
small linear-growth rate characterization of ideal crystals, in which patches are identified only 
up to translation-equivalence. This is a finer equivalence relation on patches than isometry- 
equivalence, and stronger bounds hold. The corresponding coefficient of T we obtain for 
translation-equivalence types of patches (which gives a growth rate of Nx(T) in the gap be- 
tween crystals and non-crystals) depends only on the covering radius of X, not on its packing 
radius or dimension, and we show that it is best possible in all dimensions. We conjecture that 
certain super-linear polynomial growth rates of Nx(T), while not guaranteeing that X is an 
ideal crystal, are sufficient to ensure that X has some periods — more precisely that if, for a 
sufficiently small c, Nx(T) < cT n ^^ +l for all large T then X has j independent periods — but 
no cases other than j = n have yet been proved. 
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For the repetitivity function Mx(T) (which is scale-invariant) the coefficient of T we obtain 



in Theorem f^Sj for a linear growth rate in the gap between crystals and non-crystals is the 
absolute constant |, independent of dimension. We show by examples that this cannot be 
improved beyond \. In Theorem 2.3 we also show that a coefficient slightly larger than | (but 
depending on the dimension n) in the linear growth rate of Mx(T) is enough to ensure that X 
has a non-zero period when n > 1. This coefficient depends on the packing-covering constant 



n(n) for n-dimensional Delone sets, a concept originally studied by Ryshkov [26], which is 
defined in §2. In the final section we study various properties of this constant. 

The results of this paper are analogous to results in the combinatorics of words (symbolic 
dynamics) on the lattice Z n . They can even be viewed as generalizing such results, because 
symbolic words on Z n can be encoded as Delone sets by placing points near the lattice points 
in Z n using small dispacements to distinguish the symbol types. A well known result in the 
combinatorics of words on the lattice Z — the one-dimensional case — given in Morse and Hedlund 



[21 1 is that if Ns(rn) is the number of different words of length m in a two-sided infinite sequence 
S in which A distinct symbols occur then either S is periodic and Ns(rn) < P (the length of 
the period) for all m or 

N s (m) > m + A - 1 for all m. (1.1) 

Associated to any such symbol sequence is a dynamical system with a Z-action (the closure of its 
orbit under the shift) and this result can be interpreted as a condition for this dynamical system 
to be finite. In the context of higher dimensional symbolic dynamical systems, with symbols on 
the lattice Z n and a corresponding Z n -action, the question has also been raised to what extent 
growth restrictions on the number of local symbol patterns (in rectangular patches) enforce 
periodicity, see for example Sander and Tijdeman 27, 28, [2J| and Berthe and Vuillon [^, [|] for 
a discussion of the two dimensional case. Symbolic dynamics analogues of the patch-counting 



function are called "permutation numbers" or "complexity functions" in this context, cf. [11, 18|. 
The approach of Theorem [O] can be used to show that a low enough complexity bound on the 
number of rectangular patches of such systems enforces full-dimensional periodicity. For square 
(or, more generally, cubic) patches, however, one can directly obtain an optimal result in n 
dimensions, given as Theorem |4.1| in §4; the proof is a straighforward extension of the proof of 



Morse and Hedlund 21, Theorem 7.4] of (1.1) in one dimension. 
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The analogue in symbolic dynamics of the repetitivity function is the recurrence function 



Rs(m), introduced by Morse and Hedlund [21, 22] in 1938. Given a two-sided infinite symbol 



sequence S, Rs(m) is the shortest length such that every word of length Rg(m) contains a copy 
of every word of length m that occurs in S. To aid comparison with our repetitivity function 
we shall describe results in terms of a related function Ms(m) = R,s(m) — m, which represents 
the maximum distance between the leading symbols of any two successive identical words of 



length m. Put in terms of Ms(jn), Morse and Hedlund showed [21, Theorem 7.5] that for any 



aperiodic one-dimensional repetitive sequence S in which A distinct symbols occur 

M s {m) > m + A - 1 for all m > 1 (1.2) 



and |22, p.2] that 



limsup^M>r + l, (1.3) 

m^oo Tfl 

where r = 1(1 + y/b) is the golden ratio. Moreover ( |L3|) holds with equality when S is the 
Fibonacci sequence, so represents an optimal "gap in the spectrum" result. This suggests the 
possibility that the limiting constant r + 1 in the analogous Delone set inequality on the right 
of ( |6.9[ ) in Theorem |6.2| may also be optimal. The inequality (1.2) follows from (1.1) and the 



trivial estimate Ms{m) > N,s(m) (the analogue for symbol sequences of Theorem 4.1 of @) 



yet the coefficient of m in (|1.2|) is 1, which is significantly larger than the A in our Theorem 2.2 



the analogous result for Delone sets. This seems largely due to m being restricted to integer 
values in the symbolic case. 

Some of the bounds in this paper depend on the re-dimensional Delone packing-covering 
constant ft(n), defined in §2, which was introduced and studied by Ryshkov |26j] in 1974. In 
the final section we present basic properties of this constant and the related n-dimensional 
lattice packing-covering constant kl(ti), also introduced by Ryshkov. Most of the results in 
that section are not new, but we recall them for completeness and in order to raise a question 
concerning the possible unboundedness of kl(ji) as re —* oo. This question is of interest, because 
in any dimension re with kl(ti) > 2 no lattice packing can be a densest sphere packing. 

2. Statements of results 
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Definition 2.1 A Delone set, or (r, R)-set, is a discrete set X in 1" that is uniformly discrete 
(i.e. its packing radius r by equal balls is positive) and relatively dense (i.e. its covering radius 
R by equal balls is finite). We call the values r and R the Delone constants of the set X. 

Clearly r < R, with equality only when X is a one-dimensional set of equally spaced points. 
Our definition of the Delone constant r differs slightly from other authors (see [||, [^, |26|, for 
example) who take for this value the infimum 2r of the inter-point distances in place of the 
packing radius r. 

Definition 2.2 A Delone set X in R™ is an ideal crystal if it has a full rank lattice of translation 
symmetries, i.e. X = A + F, where A is a full rank lattice in W 1 and F is a finite set. 

Note that an aperiodic set is one with no global translation symmetries, while a non- 
crystalline set may have some translation symmetries, but not a full rank set of them. 

Definition 2.3 For a Delone set X, a T-patch centered at the point x E X is 

7Mx;T) :=InB(x;T) , 

where B(x;T) is the open ball with center x and radius T, and the patch- counting function 
Nx(T) is the number (possibly infinite) of translation-inequivalent T-patches centered at points 
x of X. We use the notation "P(x; T) ~ V{y; T) to mean translation equivalence of patches, i.e. 
thatP(y;T) =P(x;T)+y-x. 

The function Nx(T) is a non-decreasing function of T and Nx(T) = 1 for T < 2r. It is 
possible that Nx(T) may be infinite for large T. 

Definition 2.4 A Delone set X with Nx(T) finite for all T is said to have finite local com- 
plexity. 

Delone sets with finite local complexity are studied in [[0^, 16, 17], where they are called 
Delone sets of finite type. The growth rate of the patch-counting function Nx(T) provides a 
quantitative measure of the complexity of X. The slowest possible growth rate for Nx(T) is to 
be eventually constant, which occurs when X is an ideal crystal, as we review in §||. Our first 
main result is that except for ideal crystals the slowest possible growth rate is at least linear 
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in T. Dolbilin et al. M, Theorem 1.3] have already shown, in the context of counting isometry 
classes of patches, that Nx(T) < c(n, r, R)T for a single value of T implies that X is a crystal, 
where 

1 

c(n, r, R) 



2(n2 + l)i?log 2 (f +2)' 

and this holds a fortiori for our function N X (T) which counts the more numerous translation 
classes of patches. But with our finer classification of patches we can obtain a larger value of 
the constant, which is independent of the dimension. 

Theorem 2.1 If a Delone set X in M. n with covering radius R has a single value of T > 
such that 



then X is an ideal crystal. 



The constant ^ here is optimal, in the sense that for any c > ^ there are Delone sets 
X £ W 1 with Nx(T) < cT, for some T > 0, that are not ideal crystals, as we show in §^. 

In [jl?], Conjecture 2.2] we put forward the following conjecture, which says that slow growth 
of Nx(T) implies that X has many independent periods: 

Period Conjecture. For each integer j = l,...,n there is a positive constant Cj(n,r,R) 
such that any Delone set X in R n with Delone constants r, R that satisfies 

N X (T) < c,-(n, r, R)T n ~ j+1 for all T > T (X) 

has j linearly independent periods. 

Weaker forms of this conjecture would be to allow Cj to depend on X or merely to assert that 
T^~ n ~ 1 Nx{T) — > as T — > oo implies that X has j linearly independent periods. Theorem [2,1| 
is a strong form of the case j = n of the conjecture. The only other case that has been proved 
to date is a strong form of the case n = 2, j = 1 in the context of doubly infinite arrays of 
symbols instead of Delone sets |l( 

The basis of the conjecture is the feeling that a pattern with less than j independent periods 
should have "(j — l)-dimensional recognizable features" at arbitrarily large scales, causing 
N X (T) to be of order at least T n ~ j+1 for some arbitrarily large values of T; but concrete evidence 
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for it is limited. We formulate it in order to focus attention on the problem, rather than to 
strongly assert its truth. The formulation above is the strongest consistent with currently known 
constraints. In particular, the conjecture cannot be strengthened to assert the conclusion when 
the inequality is satisfied for only a single value of T, as Theorem ^l] does: in [17, Theorem 2.2] 



we give examples of aperiodic sets X C W 1 , for any n > 3 and e > 0, for which there exist 
arbitrarily large radii T with N X (T) < T^ n+1 ^ 2 ^ +t . Also there is no corresponding conjecture 
(except possibly for n = 2, j = 1) generalizing the result of for the function that counts 
patches up to isometry: || gives examples of aperiodic sets X C M n , for n > 3, whose isometry 
patch-counting functions are 0(T 1+<E ) for any e > 0. 

We next consider sets that have restrictions on the distribution of patches of a given type. 

Definition 2.5 For a Delone set X, the repetitivity function Mx(T) is the least M (possibly 
infinite) such that every closed ball B of radius M contains the center of a T-patch of every 
kind that occurs in X. That is, for every T-patch V of X, B contains a point of X which is 
the center of a T-patch of X that is a translate of V . (This T-patch may extend outside the 
ball B.) 

Another way of expressing this definition is to say that Mx (T) is the largest of the covering 
radii of the sets of centers of patches of X translation-equivalent to V, taken over all T-patches 
V of X. The function Mx(T) is a non-decreasing function of T and Mx(T) = R for T < 2r. 
For large T, Mx{T) may become infinite. 

Definition 2.6 (i) A Delone set X with Mx(T) finite for all T is said to be repetitive^ 

(ii) A Delone set X is linearly repetitive if there is a constant c such that Mx{T) < cT for 
all T > 0. 

Repetitive Delone sets necessarily have finite local complexity, since the definition implies 
that Nx(T) is finite for all T. The repetitivity function Mx(T) of a Delone set provides a second 
quantitative measure of its complexity, supplementing Nx(T). Ideal crystals are repetitive and 
have M X {T) bounded. 

1 There is a parallel concept for tilings, where various terms are used: "almost periodic" in p^ ], "tilings with 
local isomorphism" in jl^, |25|, and "repetitive" in [Q. For symbolic dynamical systems the term "uniformly 
recurrent" is used. 
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Our second main result is that except for ideal crystals the slowest possible growth rate of 
Mx{T) is at least linear in T. 



Theorem 2.2 If a Delone set X in M. n has a single value ofT > such that 



M X (T) < i T, (2.2) 



then X is an ideal crystal. 



The constant | here is independent of the Delone constants r and R, which is made possible 
by the fact that the repetitivity function is scale-invariant, and it is even independent of the 
dimension n. It cannot be increased to be larger than ^ (for any n) as we show in §|6|. 



For n > 2 we can show that a slightly larger value of the constant in Theorem 2.2, while 
not necessarily guaranteeing that X is an ideal crystal, at least gives a sufficient condition for 
X to have a non-zero period. Before stating this result we need to define some constants. 

Definition 2.7 (i) For a Delone set X we call the ratio k,x = R/r the packing-covering ratio 
of X. 

(ii) The Delone packing- covering constant k(u), for dimension n, is the infimum of nx over 
all Delone sets X in M 71 . 

(iii) The lattice packing-covering constant kl{u), for dimension n, is the infimum of k\ over 
all lattices A in W l . 

These constants were intoduced by Ryshkov [Q, who actually studied the quantities \k{u) 
and ^Ki(n) (corresponding to taking 2r instead of r as the first Delone constant, as mentioned 
after Definition |2.l| ) . If we now put 

«(n) 

cm) := — r-\ 

v ; «(n) + 2 

we have: 

Theorem 2.3 If a Delone set X in W 1 has a single value ofT > such that 

M x (T)<c(n)T, (2.3) 

then X has a non-zero period. 



Clearly k(1) = 1, so c(l) = „- (reflecting the fact that in dimension 1 having a non-zero 
period is equivalent to being an ideal crystal), and K,(n) > 1 for n > 2. In §7 we present some 
results about K(n) and kl(ti), including the result of Ryshkov p^] that n(n) < 2, which implies 
that | < c(n) < | for n >2. Thus Theorem 2.3 gives at most a very slight dimension-dependent 



improvement on Theorem 2f2. In Theorem 3^ in §6 we show that one cannot increase c(n) in 



Theorem |2.3| to any value exceeding \k(tl) < 1. 

3. Bounded Patch-Counts and Ideal Crystals 

A Delone set X is an ideal crystal if and only if its patch-counting function Nx(T) is 
bounded. We have the following more precise result. 

Theorem 3.1 // a Delone set X in M. n has Nx(T) bounded with maximum value N , then X 
is a union of N cosets of a lattice in R n . 

Proof. For each T > we can classify the points of X according to their T-patches. If x and 
y have different T-patches then they have different T'-patches for every T' > T, so increasing 
T refines the classification. If Nx(T) has maximum value ./V there is a U with Nx(U) = N and 
classifying points of X according to their [/-patches gives a partition 

X = Xi U X 2 U • • • U X N 

of X such that points in the same class have identical T-patches for all T and points in different 
classes have different T-patches for all T > U. There is no loss of generality in translating 
X so that G X\. For any two points x and y we have B(0;U) C B(y;U + ||y||) and 
-B(x — y; U) C B(x;U + ||y||). When x and y are both in X\ their (U + ||y||)-patches are 
equal, so the translation by x — y that takes B(y; U + ||y||) to B(x; U + ||y||) maps the points 
of X in B(0;U) one-to-one onto the points of X in B(x — y;U). Thus x — y G X and 
Vx(x — y; U) ~ Vx(0; U). Hence x — y G X\. This shows that X\ is an additive subgroup of 
W 1 . 

Next take any xi G X x and x 2 G X 2 (if TV > 1). Then Vx{*i\ U+ ||x 2 ||) ~ Vx(0; U+ ||x 2 ||), 
and in a similar way we can use the translation by x 2 to show that xi + x 2 G X and has the 
same [/-patch as x 2 . Hence xi + x 2 G X 2 and X 2 2 X\ + x 2 . Finally, for any other point 
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Y2 G X2 we have Vx{Y2\ U + ||x2||) ~ Vx{^2', ^ + Il x 2||) and the translation by — X2 shows that 
y2 — X2 G Xi. Hence X2 -X2 C Xi so X2 = X\ + X2 is a coset of X\ in M. n . In the same way, 
each Xi is a coset of Xi. 

Since the Delone set X is the union of finitely many translates of X\ , X\ itself is a Delone 
set and, being a subgroup of M n , is a full rank lattice. □ 

Corollary 3.1 For Delone sets X, 

Nx (T) = 1 for all T <^=^ X is a translate of a lattice, and 
Nx(T) is bounded <^=^ X is an ideal crystal. 



Proof. The implications in one direction follow from Theorem 3.1. For the other direction, 
if X = A + F is an ideal crystal then the orbits of X under its group of translation symmetries 
are cosets of a lattice containing the lattice A and the number of orbits is at most the cardinal 
of the finite set F. Clearly points in the same orbit have identical T-patches for all T. □ 

4. Linear Patch-Counts and Ideal Crystals 

The object of this section is to show that if X is not an ideal crystal then there is a linear 
lower bound on the growth rate of Nx{T). The main step in doing this is the following lemma, 
which is analogous to the result in the combinatorics of words that if the word-counting function 
of an infinite sequence of symbols is the same for two consecutive word lengths then it remains 



the same for all greater word lengths (cf. our proof of Theorem 4.1 below) 



Lemma 4.1 If X is a Delone set with covering radius R and there are radii U > and 
V > U + 2R such that N X (V) = N X (U) then N X (T) = N X {U) for all T > U, and X is an 
ideal crystal. 

Proof. If x and y are points of X with Vxi^U) 9^ Vx(y',U) then clearly Vxi^V) 9^ 
Vx(y',V). So if NxiV) = Nx{U) then, for each x G X, VxfaV) is determined uniquely by 
Vx(^U). 

We next show that Vx{*-; 2V — U — 2R) is also determined uniquely by "Px(x; U). Take any 
xq G X, fixed for the moment, and consider a point c G B(xq; 2V — U — 2R). We can find a 



11 



point d G .B(xo; V — U — R) with ||d — c|| < V — R and, by the definition of R, a point x of 
X in B(d; R). Then 77(x; U) C B(x ; V) and c € 77(x, V). So 7>x(x ; V) determines P x (x; 17) 
which in turn determines Vxfal V). Since c was an arbitrary point in B(xq; 2V — U — 2R), 

B(x ; 2V - U - 2R) C |J 77(x; V) 

X 

where the union is over all points x S Vx(*-o; V — U), so the points of X in B(xq; 2V — U — 2R) 
are determined uniquely by Vx(x-o] V)- Since xo was an arbitrary point of X, 

N X {2V - U - 2R) = N X {V) = N X {U). 

Iteration of this argument successively increases V by V — U — 2R and shows that Nx(T) = 
N X (U) for all T >U. □ 



Theorem 2.1 is an almost immediate consequence of this lemma. 



Proof of Theorem 2.1. If X is not an ideal crystal then, by Corollary fO, Nx(T) is un- 



bounded. Certainly iVx(eo) > 1 for any eo > 0. Since Nx(T) is unbounded 

N X (2R + e + ei) > N x (e ) for any e x > 0, 



by Lemma 4.1, where R is the covering radius of X. Hence iV(2i? + eo + > 2. Repeating this 
argument we find N(4R + eo + t\ + £2) > 3, . . . , N(2mR + eo + ei + • • • + e m ) > m + 1 for any 
integer m > 0, where eo, • • • , e m can be chosen arbitrarily small. It follows that Nx(T) > T/2R 
for every T > 0. □ 

One can prove an analogous result in multi-dimensional symbolic dynamics, giving a con- 
dition for full periodicity of a symbol pattern on the lattice Z n . We consider symbolic words 5 
on the lattice Z n , drawn from a finite alphabet ^4, i.e. S £ A z " . 

Theorem 4.1 Let S be a symbolic word on the lattice Z n in which A different symbols occur 
and let Ns(m) count the number of different symbol patterns in S over all lattice cubes with m 
lattice points along each side. If, for some integer m > 1, 

Ns(m) < m + A- 1, (4.1) 

then S is fully periodic, with a full rank period lattice AcZ™. 
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Proof. Suppose that Ns(rriQ + 1) = Ns(mo) for some tuq. Then whenever two cubes of side 
mo with centers at x and y (where x and y are lattice points if mo is odd and half lattice points 
if mo is even) have identical symbol patterns so do the cubes of side tuq + 1 with centers at 
x+(|,|,...,2) and y + (|, 5, •••,!), since otherwise the number of symbol patterns occurring 
in cubes of side mo + 1 would be at least one greater than the number occurring in cubes of 
side mo- The same holds when (5, 21 ■ ■ ■ > 2) * s re pl &ce d by (±^, . . . , for any choice of 
signs. Consequently the symbol pattern on any cube of side mo determines uniquely the symbol 
pattern on the cube of side mo + 2 with the same center. (This cube consists of the original 
cube surrounded by an extra shell of lattice points.) Hence Ns(rriQ + 2) = Ns(m) and it follows 
by induction that the symbol pattern on any cube of side mo determines the pattern on the 
whole of Z n . Take any vector A £ Z™. Since there are only finitely many symbols there must 
be two cubes of side mo with the same symbol pattern such that one is a translate of the other 
by pX, for some Since each cube determines the same global pattern, pX is a period of 

S. Since it has periods in every lattice direction, S is fully periodic. 

It follows that if S is not fully periodic then Ns(m) is a strictly increasing function of m, 
and hence 

N s (m) > Ns(l) +m-l = A + m-1. 

□ 



Remark. The one-dimensional case of Theorem 4.1 is well-known and is due to Morse and 
Hedlund Theorem 7.4] in 1938. The proof above is an extension of the one-dimensional 
proof. The result is best possible in all cases: one-dimensional Sturmian sequences S have 
A = 2 and Ns(m) = m + 1 for all m > 1, and examples in higher dimensions are given by 
patterns that are Sturmian in one coordinate direction and constant in all other coordinate 
directions. There exist optimal examples for all other values of A too. 

5. Linear Repetitivity and Ideal Crystals 

In |L7]] we studied linearly repetitive sets, which are Delone sets X for which there is some 
constant c with Mx(T) < cT for all T. There are many examples of aperiodic linearly repetitive 



sets associated to self-similar constructions and in [17] such sets were proposed as models for 
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"perfectly ordered quasicrystals." Theorem |2,2| asserts that linear repetitivity, with a sufficiently 
small constant, is sufficient to force a set X to be an ideal crystal, and Theorem |2.3| asserts 
that linear repetitivity with a slightly larger constant, depending on the dimension, forces at 
least one period. 

Proof of Theorem |2.2| . Without loss of generality we may suppose that 0£l. Consider 
the patch V := V x (0;T) = B(0;T) n X, where T satisfies (fj), and let Ep := {y G X : 
V x (y;T) ~ V} be the set of centers of T-patches of X translation-equivalent to V. 

We first note that any vector p G Epn.B(0; |T) is a period of X. This is because x£P and 
p G Ep implies x+p G "P + p = "Px(p; T) C X. Ifx G T X (0; §T), thenx+p G B(0;T)nI = 7? 
and so, for all y G Ep, 

xGPx(y;iT) x + peP x (y;T)cI. 

since Afx(T) < \T the balls S(y; ±T) with y G Ep cover M n and therefore p is a period of X. 

It remains to show that B(0; |T) contains n linearly independent vectors pi, P2, . . . , p n of 
S-p. Let B\ be any ball of radius with on its boundary. Since jT > Mx(T), B\ contains 
a non-zero point pi of E-p. When pi, . . . ,pj have been chosen let Bi + \ be any ball of radius 
that touches the linear subspace (pi, . . . ,p,) at 0. Then there is a vector p^+i of E-p in 
Bi+i, which is necessarily linearly independent of pi,... ,Pi- Continuing in this way we find 
n linearly independent vectors pi, . . . , p n which all lie in B(0; |T), since each Bi C 5(0; §T). 
□ 

Proof of Theorem |2.3| . Let V and E-p be as in the proof of Theorem [2.2| . Since Mx(T) is 
at least as large as the covering radius of Ep the argument used in the proof of that theorem 
shows that if there are two distinct points x, y G Ep with 

||x- y|| <T-M X {T) (5.1) 

then x — y is a period of X. But the packing radius of Ep is < Mx(T)/K(n), so we can certainly 
find x, y G Ep satisfying ([O]) provided 

2M X (T)/ K (n) <T-M X (T), 

that is, provided (^^) holds. □ 
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6. Non-crystals of Low Patch Complexity 



Here we give examples that limit the extent to which the constants in Theorems 2A, [2,2| 
and |2.3| can be increased. We first point out that, for trivial reasons, the constant ^ on the 
right of ( p. 1| ) cannot be improved and the constant | on the right of (^^) cannot be increased 
beyond \ for general Delone sets. This is because for any Delone set Y with Delone constants 
r,R we have iVy(T) = 1 = ±T and My(T) = R = §T when T = 2r. When n = 1 there 
are non-crystals with ^ arbitrarily close to 1, so ^ can be made as close as we like to ^ 
and y~ as dose as we like to ^. Higher dimensional examples can be found by noting that if 
X = Y x r{E a ~ x then X is a non-crystal when Y is and has Nx(T) = Ny(T) and Mx(T) very 
close to My(T), for this T, when r\ is small. 

The growth rates of Nx(T) and Mx(T) as T — > oo are more fundamental characteristics 
of X than the sizes of these functions at a particular value of T, however, so we give below 
examples that combine near-minimal sizes of Nx(T) and Mx(T) for a particular T with small 
linear growth rates as T — > oo. In particular, these examples are linearly repetitive, which 
implies that they have a great deal of regularity in the form of the existence of uniform patch 
frequencies and of an approximately linear address map, as shown in [17, Theorems 6.1, 7.1]. 



Most of our examples are not only Delone sets but Meyer sets, where a Meyer set is a Delone 



set X whose set of differences X — X is also a Delone set, see [14|, [19]. 

We begin with one-dimensional examples and then extend them to higher dimensions by 
the technique of taking the direct product with a lattice of small mesh. 

Theorem 6.1 For each e > there exists a one- dimensional aperiodic Meyer set Y = Y(e) 
which is linearly repetitive and has the following properties. 

(i) There is a T* > with 

N Y (T*)< (l + e)^ and My (T*) < (i + e)r*. (6.1) 

(ii) There is a To > such that, for all T > Tq, 

N y (T) < (1 + e)^ and M Y (T) < (r + 1 + e)T, (6.2) 
where t = ±(1 + y/E) « 1.6180. 
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Proof, (i) As described above, to satisfy (i) it is sufficient to construct a set Y with Ky = 
-p < 1 + 2e and to take T* = 2r. We shall take for Y a set of the form 

Y = Y(a; 5) = {m + 5(ma) : m <G Z}, 

where a is irrational, \S\ < i, and (x) = x — [x\ is the fractional part of x. This set is almost 
linear, in the sense of Jl~5|| , so is a Meyer set by Theorem 5.1 of [15]. The irrationality of a 
ensures that Y is aperiodic. The packing radius ry of Y satisfies 



while the covering radius Ry satisfies 



ry > -(1 



fly < ^(1 + 1*1), 



so 



fly ^ 1 + 1*1 

Ky = < 



ry - l-\6\' 
which is < 1 + 2e when \S\ < \e. 

(ii) To bound Ny(T) for large T we note that the length of the mth interval of Y = Y(a; 5) 

is 

1 + S(((m + l)a) - (ma)) = 1 + 5a - S( [(m + l)aj - Lm-aJ ), 

so the intervals of Y~ are of two lengths and the arrangement of these lengths corresponds to 
the arrangement of the terms of the Beatty sequence (or Sturmian word) 

{[{m + l)a\ - [ma\} . (6.3) 

For investigating the symbolic dynamics of Y, the two lengths can be replaced by two abstract 
symbols, A and B say. The study of such symbolic sequences goes back to the eighteenth 
century but the emergence of a coherent body of theory dates from two papers of Morse and 



Hedlund [21, [2^] in 1938 and 1940. Coven and Hedlund ||] showed that for m > 1 there are 
exactly m + 1 distinct words of length m in such a sequence when a is irrational. Since the 
number of intervals that are contained in a given T-patch of Y, or overlap it sufficiently for the 
overlap to determine their length, is at most^] 2|_T + |<5|J, this result shows that 

N Y (T) <2[T + |«5|J + 1 < (1 + |*|)^ + 1 + 2\5\ (6.4) 



2 Perhaps more transparently, we could use instead the upper bound [T + |<5|] for th e total number of intervals 
that meet the T-patch, at the cost of slightly larger constant terms in ( |6.4[ ) and ( |6.7| ). 
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for all T > 0. 

The analogue in symbolic dynamics of the repetitivity function is the recurrence function, 
R(m), introduced by Morse and Hedlund. As in §1, we work in terms of the function Ms(m) = 
Rs( m ) ~ m i the maximum distance between the leading symbols of two successive identical 
words of length m in a symbolic sequence. It is proved in |2^, p. 2] (see also [jij, for example) 
that if {qk} is the sequence of denominators of the continued fraction convergents of a and 
Qk-i < m < Qk then M(m) = qk + Qk-i- If the continued fraction has bounded partial quotients 
(as happens for quadratic irrationals, for example) then qk + qk~i < (6 + 2)qk-\, where b is 
an upper bound for the partial quotients, so the symbolic sequence (|T^) is linearly repetitive 
and hence Y is too. To get the bound in ( |6.2j ) we take a = r, so that 6 = 1 and qk = Fk, the 
kih. Fibonacci number. Recalling that the number of identifiable intervals in any T-patch is at 
most 2 \T + 1 5 1 J, we then have 

2M Y (T)<F k + F k - 1 + \5\ (6.5) 

when 

F k - 1 <2[T + \5\\ <F k , (6.6) 

where the last term on the right of (|6.5| ) is to take account of the fact that the distance between 
two points of Y differs from the number of intervals separating them by at most \S\. Using the 
identity 

F k = TFk-i + {-) 
and the lower bound in ( |6.6[ ) we obtain 

M y (T)<(t + 1)T + ^ + (t + ^)\5\ (6.7) 

for all T > 0. If we choose \5\ < \t and T > 2(e~ 1 + e) then (0) and give fl6^). □ 

Remark. We note in passing that the set Y(t;—1/s/5) is a is a scaled version of the well 
known Fibonacci quasicrystal. It can be shown to satisfy My(T) < t 2 T + r 2 /2\/5, where the 
constant term is half the length of the long interval. In this case it can be shown that this 
bound for My(T) is optimal in the sense that there exist arbitrarily large T for which it is 
achieved. 
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It is now straightforward to extend Theorem 6.1 to the n-dimensional case by taking direct 
products with lattices of small mesh. 

Theorem 6.2 In W 1 , for each e > there exists a non- crystalline Meyer set X = X{e) which 
is linearly repetitive and has the following properties. 

(i) There is aT* > with 

^x(n< (i + e)^ and M X (T*)< (i + e )T*. (6.8) 

(ii) There is a Tq > such that, for all T > Tq, 

N X (T) < (1 + e)- and M X {T) < (r + 1 + e)T. (6.9) 
K 



Proof. We take Y = Y(e) to be the one-dimensional set constructed in Theorem 6.1 and 



X :=Y x rjZ n -\ 

where r\ is a sufficiently small positive constant. For any r], X is an n-dimensional linearly 
repetitive Meyer set which is a non-crystal, but has an (n — l)-dimensional lattice of periods. 
Because of these periods 

N X (T) = N Y (T) (6.10) 

for all T. We also have 

M x (Tf < M Y (Tf + M vZn -x{T) 2 = M Y (Tf + -^(n - 1) (6.11) 



for all T, and taking T < min(2ry,?7) in ( 6.11| ) gives 



R x <R 2 Y + jV 2 (n-l), (6.12) 



where Rx is the covering radius of X. Now ( 6.10 ), ( 6.11| ) and ( |6.12[) , together with Theorem |6.1 



(with in place of e), give ( |6.8| ) and ( |6.9[ ) when r\ is small enough, with T* and To as in 
Theorem B.lL □ 
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Remark. For the examples X(e) the constants on the right of the inequalities (|6.8| ) and ( |6.S| ) 
simultaneously approach |, 1 and r + 1. The first of these is optimal in all dimensions, by 
Theorem 531 and it is conceivable that the others are also individually optimal. We note that 
if the second value, |, is optimal then, since c(n) < \ for all n, every Delone set X satisfying 



the hypothesis of Theorem [2.3| is not only periodic but also an ideal crystal. 

We end this section by deriving attainable values for the complexity of aperiodic (rather 
than merely non-crystalline) Delone sets. Note that, since 1 < n{n) < 2 for n > 2, the constants 
m djj) below are larger than those in (i) of Theorem as is to be expected, but are also 
smaller than those in (ii) of Theorem |6.2|. 



Theorem 6.3 (i) In each dimension n > 1 and for each e > 0, there exists an aperiodic 
Delone set X = X(e) for which there is a T > such that 

N X (T) < (^K(n) + e) | and M X (T) < (i«(n) + e)r. (6.13) 

(ii) In each dimension n > 1 and for each e > 0, i/iere exists an aperiodic repetitive Meyer 
set Z = Z(e) having uniform patch frequencies and for which there is aT > such that 

N Z (T) < (l KL (n)+e^- and M Z {T) < (^i(n) + e)r. (6.14) 



Proof, (i) Take a Delone set X with kx < n{n) + e. We may suppose that X is aperiodic by 
moving one point of it by a small enough amount not to disturb this inequality but in such a 
way as to destroy any global periods. Choose T = 2rx, so that Nx(T) = 1 and Mx(T) = Rx- 
Then ( 6.13[ ) holds since 



{\^ + e ) T> (§T x + \ e ) 2rx>Rx 



(ii) Uniform patch frequencies are defined as in [17|. Take a lattice set A € M n with 
< ^i(n) + e and let 

Z = {\ + ((p(X))S : A G A}, 

where S £ M n and </> is a linear functional on R n . By choosing S small enough (for a fixed <f>) we 
can ensure that the packing and covering radii of Z remain close enough to the corresponding 
values for A for k z < Ki(n) + e to hold. Then, as in the proof of (i), ( |6.14 ) is satisfied with 
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T = 2rz- For the other properties of Z we note that, up to an invertible linear transformation, Z 
is a cut-and-project set with lattice AxZ C M"xK = M n+1 , the hyperplane x n +i = (j)(xi, . . . ,x n ) 
as "physical space" and the line joining (0, 0) to (—6, —1 — (f)(6)) as "internal space", the window 
being the half-open interval, closed at the (0, 0) end, with these points as end points. (The 
invertible linear transformation in question is projection on the K n factor, its inverse being given 
by x w (x, c/>(x).) The facts that Z is a repetitive Meyer set with uniform patch frequencies 
that is aperiodic can now be derived from Propositions 2.21 and 2.23 of |23| provided certain 
conditions are satisfied; namely, (a) the physical space contains no non-zero vector of the lattice 
A + Z, (b) the physical space is not contained in any subspace generated by n lattice vectors, 
and (c) the internal space also contains no non-zero lattice vector. Of these, (a) amounts to 
choosing <j) so that 1, 0(Ai), . . . , 4>(\ n ) are linearly independent over Q, where Ai, . . . , A n is a 
basis of A. Then (b) holds too, since the physical space, having dimension n, cannot lie in a 
subspace generated by n lattice vectors unless it contains those lattice vectors. Finally, any 
lattice vector in internal space has the form (tS, t(\ + <p{8)) with t8 G A and t(l + 4>(S)) E Z, so 
the internal space will contain no non-zero lattice vectors if S is chosen to avoid the countably 
many vectors (to — ^(A)) _1 A with A G A, to G Z and 0(A), to ^ 0. □ 

Note that while the set Z in (ii) of Theorem |6.3| is repetitive, and hence also of finite local 
complexity, the set X in (i) may have Nx{T) or Mx(T) infinite for large T. 

To give the case n = 24 as an example, using (|7.2j ) below; by Theorem 2^3 if X C M 24 
satisfies 

M X (T) < 10 ^~ 12 T « 0.4093T 
13 

for some T > then X has a non-zero period, but by Theorem [T^(ii) for any c > l/y/2 ~ 0.7071 
there exist aperiodic sets X in M 24 with uniform patch frequencies that satisfy Mx(T) < cT 
for some T > 0. 

7. Delone Packing-Covering Constant 

Recall that the Delone packing- covering constant K,(n) in n dimensions is the infimum of — 
taken over all n-dimensional Delone sets and the lattice packing-covering constant kl(u) is the 
infimum of y over all n-dimensional lattices. In this section we prove various facts about n(n) 
and kl(ti). We first show that these infimums are attained. 
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Theorem 7.1 In each dimension n there exists a Delone set X whose Delone constants (r, R) 
satisfy = n{n), and a lattice A whose Delone constants satisfy -p = k,l(ji). 

Proof. For each m G N let 

5 m = {ScR R : |x| < m and |x - y| > 1 for all x, y G S}. 
Then each S G S m contains only finitely many points and the Hausdorff metric 

d(Si, S2) = max min |x — y| + max min |x — y| 

xe5i yeS 2 xg5 2 yeSi 

makes S m a compact metric space. After suitable scalings we can find a sequence {X{\ of 
Delone sets in M. n having r > 1 for all i and R tending to n(n) as i — > 00. The sets Aj n -B(0; 1) 
belong to the compact space Si, so we can find a convergent subsequence 



which we can choose so that d(Su, Sn) < 1 for all i. Similarly the sets Xu n -6(0; 2) belong to 
S2, and we can find a subsequence 



with d(S2i,S2i) < 1/2 for all i. Continuing in this way, for each j G N + we can find a 
subsequence {Xji}^ =1 of {Xj—i > i}f2. 1 such that 



and R < n{n), and by the definition of n(n) there is equality in both places. 

The proof for the lattice case is similar, noting that a limit of lattices that is a Delone set 
is a lattice. □ 

We next give upper and lower bounds for «(n), which are due to Ryshkov [pt]]. 
Proposition 7.1 In each dimension n the Delone packing-covering constant «;(n) satisfies 



{Su} = {XunB(0;l)} 



{S 2i } = {X 2i n B(0; 2)} 



d{Xji n B(0;j),Xji n B(0; j)) < - for all i G N + . 
The sequence of Delone sets {Xji}JL 1 now converges to a unique set X C 1" that has r > 1 




(7.1) 
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Proof. The lower bound for n{n) follows by combining Rogers's upper bound for the packing 
density Theorem 7.1] and Coxeter, Few and Rogers's lower bound for the covering thickness 



[33, Theorem 3.4] of equal balls in W 1 . The former asserts that for every e > the number 
of balls of radius r that can be packed into a ball of radius I is < (1 + e)a n (l/r) n when I is 
large enough, where a n is the proportion of an re-dimensional regular simplex that is covered 
by balls centered at its vertices and reaching to the mid-points of its sides. The latter asserts 
that the number of balls of radius R it takes to cover a ball of radius I is > (1 — e)r n (l / R) n 
when I is large enough, where r n is the proportion of the simplex (counted with multiplicity) 
that is covered by balls centered at its vertices and reaching to its centroid (these latter balls 
just cover the simplex). It follows from these bounds that the packing radius r and covering 
radius R of any set ICR™ satisfy 

l/n , 

2re 




re + 1 

the expression on the right being the ratio of the circumradius to half the edge length of a 
regular n-simplex. 

To establish the upper bound for «(n), take a packing or W 1 with unit spheres, to which no 
sphere can be added without overlap. The sphere centers then form a Delone set X with r > 1. 
We claim that R < 2. If not, and R > 2, then there would exist some point at a distance at 
least 2 from all sphere centers, and a unit sphere placed at this point would overlap no sphere 
in the packing, a contradiction. □ 

The lower bound in Proposition [O] is sufficient to establish that k(1) = k_l(1) = 1 and 
k(2) = kl(2) = 2/\/3, with the optimal configurations in dimensions 1 and 2 being Z and the 
hexagonal lattice A2. This gives c\ = | and C2 = |(\/3— 1) ~ 0.3660. In three dimensions we 
have the bounds ^ 

i/| « 1.2247 < k(3) < k l (3) = y| « 1.2910, 

the upper bound coming from the body-centered cubic lattice A%. In twenty- four dimensions 
we have the remarkably close bounds 

~ 1.3856 < k(24) < k l (24) < \[2 1.4142, (7.2) 

with the upper bound coming from the Leech lattice. (The values of the packing density and 
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covering thickness of the Leech lattice are given in 0].) Note that n{n) and K^(n) are probably 
not monotonically increasing in n. 



The lower bound in Proposition [7JJ also yields the asymptotic bound 

liminf n{n) > V2. 

n^oo 

This asymptotic bound can be further improved by using the current best upper bound for 
sphere packing density (the Kabatiansky-Levenshtein bound) as n — > oo, instead of the Rogers 
bound, which yields 

liminf K(n) > 2 a5990 1.5146, 

n— >oo 

cf. @, Chapters 1 and 9], @, Theorem 8.2]. 

The upper bound of Proposition |7[l| shows that the constant c(n) = K ^"_| 2 — 3> nence 
Theorem |2.3| gives at most a very slight dimension-dependent improvement on Theorem |2.2|. 



The lattice packing-covering constant kl(u) satisfies the lower bound in (7.1) a fortiori. In 



[26 1 Ryshkov gives the upper bound 



I \ <r 71 + 2 

(Recall that Kx(n) is, by definition, twice Ryshkov's constant.) This is known to hold with 
equality in dimensions 1, 2 and 3. There seems to be very little known about upper bounds for 
Ki(n), and we wish to raise the following question.^ 

Question. Is there a constant C such that kl{ti) < C for all n, or is Ki(n) unbounded as 
n — > oo ? 

One reason for interest in the question is that in any dimension n for which ni{n) > 2 (should 
such dimensions exist) the densest sphere packing with equal spheres cannot be a lattice packing. 

Acknowledgment. We are indebted to V. Berthe for comments and references concerning 
Theorem |6.1| . 
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